Charges on Strange Quark Nuggets in Space 
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Since Witten's seminal 1984 paper on the subject, searches for evidence of strange quark nuggets 
(SQNs) have proven unsuccessful. In the absence of experimental evidence ruling out SQNs, the 
validity of theories introducing mechanisms that increase their stability should continue to be tested. 
To stimulate electromagnetic SQN searches, particularly space searches, we estimate the net charge 
that would develop on an SQN in space exposed to various radiation baths (and showers) capable of 
liberating the SQN's less strongly bound electrons, taking into account recombination with ambient 
electrons. We consider, in particular, the cosmic microwave background , radiation from the sun, 
and diffuse galactic and extragalactic ultraviolet backgrounds. The largest charge, for the settings 
considered, develops on a solar system SQN exposed to a solar X-ray flare. A possible dramatic 
signal of SQNs in explosive astrophysical events is noted. 



I. INTRODUCTION 

A generation has passed since Witten [l[ suggested 
that strange quark matter (SQM) might be the ground 
state of all familiar matter and de Rujula and Glashow 
@ set down a list of methods to search for strange quark 
nuggets (SQNs) in various mass ranges. There have been 
significant efforts using most of those methods, but no 
SQNs have been found. An area that has received, we 
believe, little attention is the possibility of exploiting the 
charge on SQNs in space that should be caused by ra- 
diation baths liberating electrons from the SQNs. How 
challenging such exploration would be depends on how 
large such charges might be. We estimate those charges 
under various conditions, including those at the time of 
cosmic recombination, those in today's cosmic microwave 
background (CMB), those in the solar neighborhood from 
the quiet sun and from an X-ray flare, and those in the 
diffuse galactic and extragalactic radiation backgrounds. 

We begin, in this section, with a brief review of SQM 
and SQN basics. In Sec. II, we give the equations on 
which our numerical estimates are based and the approx- 
imations that go into deriving them. We make simple, 
conservative approximations in estimating electron bind- 
ing energies, wave functions and cross sections for large 
and small SQN masses and radii. In a separate pub- 
lication, we will show, in more detail than needed for 
the first estimates here, the wave functions in the tran- 
sition between these two regions 0. In Sec. Ill we give 
the numerical results. In Sec. IV, we briefly discuss the 
results and their implications for the feasibility of space- 
based or space-directed electromagnetic SQN searches. 
We conclude, in Sec. V, by collecting evidence, pro and 
con, on SQN existence to provide the reader a context in 
which to decide whether to devote time to exploiting our 
results for the charges on SQNs in space. 

In 1984, Witten [l[ (and Bodmer Q earlier) considered 
systems of up, down and "strange" quarks, pointing out 



that they would have the same attractive potential en- 
ergy as systems of just up and down because the force be- 
tween two quarks does not depend on their flavor. They 
would, however, have about 10 percent less kinetic energy 
because the Pauli exclusion principle would not force the 
quarks into as high kinetic energy states as in the case 
with just two kinds of quarks. 

Soon after this seminal work suggesting that SQM 
might be the lowest energy state of familiar, baryonic 
matter, Farhi and Jaffe [f| worked out the basic nuclear 
physics of strange quark matter (SQM) within the MIT 
bag model [||, and de Rujula & Glashow Q identified sev- 
eral observations that might lead to discovery of SQNs. 
The latter proceeded from the basic relation for energy 
loss 



dE/dt = -Trr%p M v^f 



(1) 



where vjy is the speed of the (spherical) nugget, rjv its 
radius, and pu the density of the material through which 
it is passing. Equation (1) just says that the SQN must 
lose the energy needed to make the stuff in its way move 
as fast as it, the SQN, is moving. They estimated SQN 
mass ranges to which various sensors might be sensitive, 
for example the Earth network of seismometers being sen- 
sitive to masses over a ton or telescopes being sensitive 
to light emitted by SQNs entering the atmosphere. They 
estimated the upper detectability boundary of mass re- 
gions by computing the mass M at which events would 
become too rare for the detector system if the galactic 
dark matter (DM) density, pdm, were all in the form of 
SQNs of mass M. Roughly, we have, for events per unit 
time, 



dN/dt = v N (p DM /M)irr 2 d 



(2) 



where vn is again the speed of the SQN, the dark matter 
density in the solar neighborhood is pdm = 5 x 10~ 25 
gem - and r c i is the radius of the detector system 
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through which it is passing. Equation (2) gives minimum 
detectable mass or speed in terms of the dark matter den- 
sity limit on the abundance of SQNs of some one single 
mass. 

An important recent development is work by Alford et 
al. 01 showing that, for SQM in bulk, Cooper pairing, 
the basic phenomenon of superconductivity should occur 
with quarks. The pairing could take different forms. For 
a useful review, see Alford et al. [U Most likely at high 
density, perhaps, would be pairs of quarks with equal 
and opposite momentum, antisymmetric in spin, flavor 
and color: color- flavor locked (CFL) pairing. The quark 
matter of such an SQN would be electrically neutral in its 
interior, but would have net positive charge on its surface, 
with the total (quark) surface charge proportional to the 
area. It would have quark charge Zq = 0.3A 2 / 3 (where 
A is one third the number of quarks) for an SQN with 
m s = 150 MeV (see [§[), balanced by electrons feeling 
a potential that is a step function at the SQN radial 
boundary and Coulomb beyond it. We will discuss this 
model further below. 

There are two potential sources for SQM. It might 
have been produced primordially, in particular in a 
phase transition in the early stages of the big bang. 
That was Witten's original thought, but it is, at best, 
controversial. The issue is that an SQN formed at high 
temperature needs to cool. If it does that by evapora- 
tion, the SQN disappears. Alternatively, it might cool 
by neutrino emission in which case it survives. There 
are experts and arguments on both sides (see [l(J for a 
review and Section II. D below for a brief discussion). 
Primordial SQNs would likely all be very close in mass 
assuming limited colliding after formation. The second 
potential source is "neutron stars (NS)." If SQM is 
the lowest energy state of matter, it is expected that 
Type II supernovae would likely rise to high enough 
temperature to cool into SQM. If one did not, local 
quantum fluctuations would likely soon cause a global 
transition of the NS to an SQS (strange quark star). 
Binary SQS systems would in time spin down, collide, 
and the galaxy would gain a population of SQNs of 
varying masses from the fragments. 

There is significant, ongoing activity in searching for 
SQNs in various mass ranges. See Finch [ll[ for a re- 
cent, brief review. One important effort is the Alpha 
Mass Spectrometer experiment [l2j which was due to be 
flown by the space shuttle to the space station in a cou- 
ple of years. It would have been able to detect, and 
to distinguish from cosmic rays, light SQNs that would 
not penetrate the atmosphere. However, the Columbia 
accident and the need to retire the shuttle fleet and to 
complete the space station have prevented AMS launch 
as scheduled [13[. A second is the ICE CUBE neutrino 
detector [l4| being installed in Antarctica. It will have 
photomultiplier tubes to detect the products from colli- 
sions of (weakly interacting) neutrinos with the electrons 
and nuclei of the ice. It will also have acoustic detectors 



TABLE I: Some Strange Quark Nugget Searches. 
Experiment/Observation Mass Range (g) Result 



AMS" 
RHIC 
Mica Tracks 6 
ICE CUBE C 
Seismometers: 

Future Lunar d 

Apollo 6 

USGS Reports c 



10 



-24 _ 

< 3 x 1(T 21 

-20 



1(T 2() - 10" 
10" 3 - 10~ 2 



W* - 10 b 
10 4 - 10 6 
10 6 - 10 s 
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not done 
not found 
<< Pdm 
not done 

not done 

< Pdm/W 

< Pdm 



a Sandweiss[]J. 
b Price [l5l| 
c Spiering[l4j. 
d Banerdt et ai.llil . 
e Herrin et aUl9H. 



to be able to identify tracks made by SQNs. It should 
be able to detect SQN masses up to as much as about a 
gram. This is the highest mass for which Eq. @ gives a 
dozen or more events per year in a kilometer-sized region. 

Past searches have included examination of tracks in 
mica by Price at Berkeley [l5[ . They have also included 
two cases of NASA work, in 2002, with evidence that 
two neutron stars were actually strange quark stars [IB]. 
However, it was later concluded that alternate explana- 
tions for the observations were more likely [l7| . 

Looking for evidence of SQNs was a prime objec- 
tive of the Relativistic Heavy Ion Collider (RHIC) at 
Brookhavcn National Laboratory. Since RHIC collides 
gold nuclei with gold nuclei, SQN masses up to two gold 
nuclei could, in principle, be produced. The experiment 
found none [121 ]. However, the binding energy per quark 
of SQM would increase with increasing numbers of quarks 
Q , so it would not be surprising if systems with hundreds 
or even many thousands of quarks were not found, but 
larger assemblies were found. 

Selected searches are summarized in Table I. In the ta- 
ble, a few major searches are listed (we modestly include 
ours [3 [H|])j along with the mass range to which they 
will be/are/were sensitive and the result, where there is 
one, in terms of the inferred SQN density in our region of 
the galaxy. An interesting SQM space search using the 
equipment being deployed to monitor near Earth aster- 
oids has been proposed by Horvath [20j. Our methods 
below should apply to other, similar models such as that 
of Zhitnitsky [2l| on cold dark matter as compact com- 
posite objects. 

Our ultimate goal is electromagnetic space tests, in 
several different settings, designed to discover SQNs or 
to falsify current SQN models. The near-term goal of 
this paper is to determine the charges that would develop 
on SQNs in these settings. Each setting is characterized 
by a photon distribution and an ambient electron dis- 
tribution. In the next section we discuss the settings, 
models of SQN structure, and the formulary. In the fol- 
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TABLE II: Settings. The quantity zm is the red shift at 
formation of the first stars; vsmbh is the radius of the super- 
massive black hole at the Milky Way center; (I)GDR is the 
(inter) galactic diffuse radiation and rg is the solar radius. To 
is the temperature of the quiet sun and T x of solar X-ray flux. 



Location 


Radiation Source 
Intergalactic 


Galactic 


Solar 


Intergalactic 


CMB: (1 + z)T 








+IGDR 






Galactic 


IGDR 


GDR 






ZlII > z 


r > r S MBH 




Solar 


IGDR 


GDR 


r > rs 








T or T x 



lowing Section (III) we present our results. In Section 
IV, we day dream of possible electromagnetic detection 
schemes. In the final Section (V), we give a brief sum- 
mary of evidence for and against SQM existence. 



II. THE FORMULARY AND ITS 
APPROXIMATIONS 

A. Settings 

We consider extragalactic, galactic and solar system 
SQNs. In each case, we specify photon and electron en- 
ergy and angle (i.e., radiation bath - isotropic, or unidi- 
rectional shower) distributions. We go back in red shift 
z to recombination at z=1089.. 

In Table II we summarize these settings. r sc is the dis- 
tance from the sun to the center of the galaxy, vsmbh is 
the radius of the supermassive black hole there, and r$ is 
the solar radius. For the cosmic microwave background 
(CMB) wc take just blackbody radiation at temperature 
(1 + z)Tq = T where To = 2.75K, today's temperature. 
Galactic diffuse radiation (GDR) and today's X-ray in- 
tergalatic diffuse b ackg round radiation (IDBR) are both 
given by, inter alia, [22f | and [23[ . Let the ionization of the 
SQN be Zm- For each of the settings, we determine Zjy 
by setting the rate of ionization Z+ by the photon bath 
equal to the rate of electron capture Z_ from the ambi- 
ent electron bath. For each setting we need estimates for 
the speed and density of free electrons (or even hydrogen 
atoms in cases where Zn is large enough to rip off an 
atomic electron) as well as for the energy distribution of 
the photon bath or shower. 

SQN settings considered are the following. (1) For in- 
tergalactic SQNs we have three settings: the CMB with 
z = 1089 and z = as well as the IGDR with z = 
along with the dominant electron bath. (2) There are 
three settings for SQNs located in the solar system with 
charges determined by the competition between the pho- 



tons of diffuse radiation, sunlight or X-ray flares and the 
electron shower of the solar wind. (3) Finally, we also 
have SQNs at the center of the galaxy (COG) subject to 
the local radiation and electron baths there. 



B. SQN Structure Models 

We consider cases with SQN structure such that the 
quarks are free to move (e.g. see Q)- Then, quarks and 
electrons form Debye clouds which migrate to the surface, 
minimizing the energy of the system and leaving behind 
an electrically neutral system of quarks. We make three 
major approximations in computing the binding energies 
of the electrons. First, we approximate the potential en- 
ergy of the least bound electron, in the case that Zm 
electrons have already been ionized, as in Fig. 1, that is 

V z {r) = -Z N e 2 /r, r > r N ; 

= oo, r < tm- (3) 

Second, where convenient, we replace the infinite bar- 
rier of Eq. ([3]) and Fig. 1 by the potential of Fig. 2, an 
apparently violent approximation that doesn't actually 
make much difference at all except perhaps in the tran- 
sition region, around nanogram masses, from atom-like 
SQNs (SQN radius less than Bohr radius) to large SQNs 
(opposite inequality) as discussed below. The work of 
reference 0] is designed to cover that transition region. 
Third, we consider only s-wave electrons moving in the 
effective field of all their predecessors and the quark lat- 
tice. This is a conservative approximation: by omitting 
angular momentum we are underestimating kinetic en- 
ergy and hence underestimating the number of electrons 
that are ionized in the given setting. 

These approximations are conservative with respect to 
our goal of identifying ways of detecting charged SQNs 
in space. All three overestimate the amount of binding 
by omitting the centrifugal barrier, the electron-electron 
repulsion and the Pauli principle for the electrons. There- 
fore, they underestimate the amount of charge the SQN 
will acquire, and thus underestimate its detectability. 

The result of the three approximations is to change the 
potential to that of Fig. 2. The calculations below are 
based on that potential. 

An important case is color-flavor locked (CFL) pairing 
0, in which the u, d, and s quarks pair symmetrically. 
Then, everywhere in the bulk there are the same number 
of quarks of each kind and hence charge neutrality of the 
bulk of the quark lattice. In this case the charge lies 
at the surface where the longer Compton wavelength of 
the (low mass) u and d quarks requires that the edge 
quarks be u and d unbalanced by any of the much more 
massive s quarks. This gives a net total quark charge on 
the order of Zt ~ (M/m p ) 2 / 3 , where m p is the mass of 
the proton. The full system has, in vacuo, in addition 
enough electrons so that the total charge of an isolated 



4 



V(r) 



e z Z N 



2r B 3r N 4ra ta, £j=„ 



FIG. 1: Potential for least bound electron. 
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FIG. 2: Approximation to potential for least bound electron. 



SQN is zero, and again the potential is that of Fig. 1 
which can be approximated here by that of Fig. 2. The 
discussion in this paper is meant to apply directly to the 
CFL model. 

In Fig. 1, the cramped quarters might appear, by the 
uncertainty principle, to make the kinetic energy appre- 
ciable, therefore yield much less binding than the po- 
tential of Fig. 2, and hence give higher Zn values and 
easier detection. However, the quarters are not really 
so cramped: the distance from r = r/v to r such that 
V(r) = y(rjv)/2 is rjv- This means that, for a given 
setting, the kinetic energy must fall with increasing SQN 
mass as M -2 / 3 while the potential energy falls only as 
M -1 ' 3 . This occurs for the high M region, rjv > o>b/Zn, 
with as the Bohr radius for hydrogen. In that case, 
we show below, Eqs.© and/or (|10[) (or, numerically, in 
Figure 3), that Zn ~ b ~ M 1 / 3 . If the kinetic energy 
becomes negligible, halving it doesn't much matter. 

The opposite case, low- mass atomic model, is < 
cib/Zn- There, the SQN is sufficiently small that the 
least bound electron is atomic in nature, in the sense 
that the overlap of its wave function with the SQN is 
small. Thus, the effect of the third approximation is 
small. We can get a sense of the effects of the other 
approximations by considering the variation in binding 
energy in the periodic table. Our approximation treats 
all the SQNs (with rjy < cib/Zn) like hydrogen. In real 
life, the most tightly bound last electron in the ground 
state of any atom in the periodic table is that of helium 
with binding energy about twice that of hydrogen while 



the least tightly bound, ground state electron is that of 
cesium with binding about one third that of hydrogen. 

There are multiple possibilities for SQN structure ac- 
cording to its surface tension (See Alford [24j). For sur- 
face tension below a critical value, large CFL SQNs de- 
cay into smaller ones while for surface tensions in the 
neighborhood of the critical value, an SQN develops a 
crystalline crust with smaller SQNs of sizes around a De- 
bye length immersed in a gas of electrons. We assume 
below that neither of these two cases obtain. Rather, we 
assume: there are only CFL SQNs with surface tension 
sufficient to ensure a neutral baryonic bulk; and there is 
a uniform surface of u and d quarks retaining sufficient 
electrons to give total charge zero in the absence of ex- 
ternal ionizing effects. We then estimate, for different 
external settings, the charge that would develop on such 
an SQN. The upshot is that our results cover SQNs of 
masses for which SQNs may not exist or may not be sta- 
ble. Even if SQM is the lowest energy state of (baryonic) 
matter and SQNs do indeed exist - they are in other 
mass ranges. Recent work on SQN surface tension has 
been done by Oertel and Urban [25[ 



C. The Formulary 

We equate the rate, Z + at which electrons are ionized 
by photons in the radiation bath (or unidirectional "ra- 
diation shower" in some cases ) to the rate at which they 
are replaced by capture from the ambient electron bath 
or shower, Z_. Our expression for Z + is 



Z+ = Tib 2 J dE N 7 (E) 

J Z N e 2 /r N 

N e {E B < E) cr( 7 + SQN -> e + SQN), l] (4) 

where: b =max [r/v, o-b/Zn\) o-b is the conventional 
Bohr radius [?ic/(amc 2 )]; r/v is (still) the SQN radius 
[(3M/47T/9) 1 / 3 = rjy]; N 1 {E)dE is the flux density of 
photons with energy between E and E + dE; N e is the 
number of electrons per unit area with binding energy 
E or less; a is the particle physics cross section for ion- 
ization and can be evaluated using standard relativistic 
quantum mechanics. Here and below, Eb = Z^e 2 /b is 
the binding energy of the least tightly bound electron 
in the SQN. The square bracket under the integral sign 
is the probability of the photon, when it hits the SQN, 
actually liberating an electron. Because a probability is 
needed, when that product is greater than one, it must 
be replaced with one. 
Similarly, wc have 



N 



,-E E 



v e (E)n e (E) l + f e (E,Z N ) 



h(E)g(e + SQN -> SQN + X,E)dE (5) 
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where h(E) is the distribution of incoming electron en- 
ergies and g is the probability that the incoming elec- 
tron will be slowed and captured times the effective area 
(of the possibly macroscopic SQN) over which capture 
is possible. We integrate over the distribution of incom- 
ing electron energies; v e (E) is electron speed; n e (E)dE 
is the number of electrons per unit volume with energies 
between E and E + dE. The function f e is given by 

f e (E) = ahcZ N /(r N E) (6) 

fe{E) is the enhancement of the effective cross sectional 
area of the SQN due to the fact that the SQN charge 
focuses the incoming electrons. It is the electrostatic 
analogue of the gravitational focusing factor discussed in 
detail on page 541 of Binney and Tremaine [26|. It was 
introduced by V. S. Safronov [2tJ as an enhancement of 
the accretion rate of planetesimals in a dusty disk. Note 
that the coefficient on the right-hand side in Eq. ([J} is 
irb 2 while the one in Eq. ([5]) is irr 2 N , (with b the larger of 
rjy and othc / (E b Z n)) ■ A photon can eject an electron 
out to b while an electron must penetrate to to be 
captured in our approximation. 

Scattering off the bound electrons or the much more 
numerous quarks and the absence of rclativistic electrons 
at the Z values under consideration make it a reasonable 
approximation to assume that all electrons up to total 
energies significantly higher than m e + Eb are captured. 
In the spirit of this approximation we set the product, 
hg, equal to 8{E — E e ), where E e is the average ambi- 
ent electron energy, and replaced the whole integral by 
n e v e (l + f e ). This approximation would need to be re- 
visited were we to consider an explosive event such as a 
supernova with MeV energies. 

We solve the equation Z + = Z- in the various settings 
described, letting M range from 10~ 21 g to 10 +30 g. We 
continue to make the conservative approximation that 
the kinetic energy is negligible compared to the potential 
energy for large and hence use, for the binding energy, 
Eb ~ Z N e 2 /r N for r N > a B /Z N . 

For large Zjv, each SQN will be surrounded by an ex- 
cess of unbound electrons and there will be screening that 
will affect both electron capture, Eq. ([5]), by limiting the 
distance over which there is an attractive force, and elec- 
tron liberation, Eq. (fj|, by impeding escape to infinity. 
In the first approximation of this work, we do not at- 
tempt to estimate the size of these effects, but note that, 
because they have opposite effects, there is a possibility 
of first order cancellation. 

We have solved the equation Z + = Z_ for some values 
of the parameters, and with some (further) approxima- 
tions. The most important of the approximations is re- 
placing the product N e a in Eq. ^ with unity. This is a 
good approximation because, with Zt ~ M 2 / 3 , we have 
N e = Z T /(4nb 2 ) ~ few x 10 25 cm -2 — independent of 
mass M. Since a ~ 10~ 20 cm 2 , the probability of a pho- 
ton liberating an electron from an SQN is of the order 
of one if the photon's impact parameter lies within the 



effective cross section from the center of the nugget, and 
if the photon is sufficiently energetic. This approxima- 
tion, nb 2 for the capture cross section, should be valid for 
electron energies (non-relativistic) of the settings in this 
paper (see Table III). 

The final result of these approximations in Equa- 
tions ([U and ([SI) is that Z + = Z- reduces to 

nb 2 cF y (E > E B ) = irr 2 N n e v e (l + f e ) (7) 

where: b is given just below Eq. Q; F^(E > Eb) is the 
number of 7 's per unit volume with energies greater than 
Eb] Eb is the binding energy of the least bound remain- 
ing electron in the SQN; n e is the density of ambient 
electrons; v e is their (average) speed; f e is the (classi- 
cal) focusing factor of Eq. ((HJljand barred quantities here 
and below are ensemble averages. For four settings - 
z=1089, z=0, quiet sun, and X-ray flare — we use a 
thermal distribution. For three — intergalactic, near the 
sun, and milkyway center diffuse radiation backgrounds 
(DRB) — we first used a non-thermal, power law spec- 
trum approximating the X-ray spectrum graph [22j ■ See 
also [23[. However, in the end it was the UV spectrum 
which was important. The values are from page 143 of 
Kolb and Turner [28| . The number of photons per unit 
volume with energy E 1 = E > Eb, denoted F(E) is given 
by 

F(E) = -c- 1 (p + l)- 1 AEV + V (8) 

with j3 = —3.05; Log(A) = 7.1 for intergalactic diffuse 
radiation spectrum. 

Our approach to the SQN clcctrosphere is complemen- 
tary to the classic treatment of Alcock, Farhi and Olinto 
[291 ]. They consider the electrosphere as a continuous 
medium in the Thomas- Fermi approximation [3(| ■ We 
consider the electron with the least binding energy in 
the mean field (of the positive lattice plus the remain- 
ing electrons) approximation. Additionally, we do not 
make use of the fascinating work of Usov [31| and oth- 
ers [HI on SQN cooling by thermal pair creation based 
on Schwingers work on pair production in strong elec- 
tromagnetic fields [HI]. We assume that the SQNs are 
cooled well below the temperatures at which thermal pair 
production becomes important. Our consideration below 
of pair production is for a cold, even if highly charged, 
SQNs. 

III. RESULTS 

Figure 3 gives results for the equilibrium values of Zn 
for seven selected settings. The seven curves, proceeding 
from top to bottom, are 

1. the sun shining on an SQN near the Earth during an 
X-ray flare; 
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Log(Z N ) 




Log (M) 



FIG. 3: SQN charge Z N (M). From the top, the 7 curves are 
in the order described in the text at the beginning of Section 
III, and used for the 7 rows of Table III. Mass, M, is in grams. 



2. the diffuse background radiation (DBR) shining on 
an SQN in the intergalactic medium (IGM); 

3. the (quiet) sun shining on an SQN in solar orbit near 
the Earth; 

4. an SQN in the primordial universe at recombination; 

5. the galactic diffuse radiation (GDR) at the center of 
the galaxy (COG) shining on an SQN located near the 
center; 

6. the GDR shining on a solar system SQN near the 
Earth; and 

7. an SQN in the CMB today (ignoring the DRB and 
any other radiation). 

In the Log-Log plot of Fig. 3, the third, fourth, and 
fifth curves from the top are close to each other. The 
little, relatively flat tails on the left in Fig. 3. give a 
rough approximation to the behavior in the mass region 
in which there is a transition between the "atomic model" 
(rjy < as/Zjv) and the M 1 / 3 behavior for > <ib/Zn. 
In Fig. 3 we have truncated curves where the computer 
program gives Zm{M) < 1 

Table III gives the parameters for the radiation and 
electron baths and showers in the same order as just de- 
scribed. In two cases where parameters arc unknown, 
we use our best estimates. One is the electrons in the 
IGM. We know that there about 10~ 6 baryons per cm 3 
there. Recent sightings of "warm-hot intergalactic mat- 



ter" - WHIM (See, for example, Paerels et al. [34| for 
a brief review and further references) indicate potential 
abundance of WHIM in the IGM as high as fifty per- 
cent. We assume ten percent in a rare burst of caution. 
The energy distribution of the electrons in the WHIM is 
not known. For the calculations, we take a generic 100 K 
for it. The second case where we estimate parameters is 
that of the center of the galaxy (COG), where conditions 
are not well known outside the immediate neighborhood 
of the supermassive black hole (SMBH). We define the 
COG to mean the first 100 pc. This gives a volume large 
enough to be expected to have a significant number of 
SQNs, if they exist, as well as a number of baryons or- 



TABLE III: Parameters for the seven curves of Fig. 3. E is 
in eV, n e has units of cm -3 , v e has units of cm/sec. F(E) 
is the function defined by Eq. ((8} for the non-thermal distri- 
bution cases. F(E) = AE~ 3 ' 05 ~f /(cm 2 seV) where average 
temperatures are not known, we assume 100 K. 



SQN Location 


Radiation 




n e 


v e /W 5 


Solar X-ray 


T = 10 3 eV 




7 


500 


Flare at 1 AU a 










IGM 6 


F(E); Log(A)= 


=7.1 


HP 7 


60 


Near Quiet Sun c 


T = 0.5 eV 




7 


500 


IGM Pre Recombo 


1 CMB T = 0.26 eV 


287 


300 


Galaxy Center e 


F(E); Log(A)= 


=11.57 40 


60 


Near Sun(DRB) / 


F(E); Log(A)= 


=7.64 


7 


500 


Today (CMB) 9 


T = 2.75K 




lO" 7 


60 



"Ref. [36j p. 423; p. 409; p. 409. Here and below, 1st entry is 
radiation; 2nd electron density; 3rd electron speed. 

6 Ref. [23 p. 143; Ref. [ll pp. 662-3, assume 0.1 ionized H; 
assume T = 100 K. 

p. 396; p. 409; p. 409. 
same for electrons. 

p. 570; Ref. [36| p. 571 extrapolated to 100 pc; assume 




/Ref, 
SRef. 



p. 143jRef. [3J p. 409; p. 409. 
Ref. [H pp. 662-3 ; pp. 662-3. 



ders of magnitude larger than the number one would have 
from the SMBH's 10 6 solar masses, but small enough to 
have approximately similar conditions throughout. As- 
suming r~ 2 behavior, there are 10 9 solar masses of pro- 
tons in the COG, (see [35[) , and assuming 10 percent of 
that is ionized, gives us a value to use for the electron 
density. Again, we take a generic 100K for the electron 
temperature. 

Based on Table III, one can understand the order of 
the curves in Fig. 3. The largest SQN charge, among 
the settings considered, is in the X-ray shower from a 
solar X-ray flare. The neighborhoods of supernovae and 
other explosive astrophysical events, however, might well 
be even better, depending on the extent to which assets 
like the Swift satellite are able to get sufficient data, in 
seeing SQN effects. The CMB by itself is, today, the 
worst, but it is perhaps interesting that, for sufficiently 
massive SQNs, it would contribute. 

In Fig. 4, we give the times (in years), T eq (M), neces- 
sary to reach Zjy [r eq (M) = Z^/Z + = Z^ / ZJ\ for (just) 
the dominant contribution in each location. The results 
for Zm do not vary with distance from the sun because 
both the radiation and electron showers fall off like r~ 2 . 
However, the times to reach Zn do vary with r like r 2 . 
The times are about a day and a half for a nanogram, 
and fall as M -1 ' 3 as shown in Fig. 4. See the discussion 
of Equation (9) below for a derivation of this result. 

Roughly, there are three thick lines in Fig. 4. The 
middle line, with the intermediate times is composed of 
three individual lines hard to distinguish in a log-log plot. 
They are the settings 1, 3, and 5 discussed above and arc 
the three solar system settings (near the Earth). The 
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Log(r eq (M) ) 

10r 




Log (M) 



FIG. 4: Time in years to reach equilibrium Zn(M). Top 
double line is IGM and CMB today; next, triple, line is the 
solar system at Earth's distance from sun during X-ray flare, 
from GDR, and for quiet sun. The bottom double line is 
galactic center and universe at recombination. 



TABLE IV: Times and Binding Energies. The settings are 
grouped into three sections: the three solar system ones to- 
gether, the two intergalactic ones together and the universe 
just before recombination which is close for no apparent rea- 
son to COG. These settings are not the same as in Fig. 3 and 
4 or Table III. 



Setting 




(y) 


E B {eV) 


E B (eV) 








Al > 10- 10 g 


i(r 21 g 


Galactic Center 


38 x 10" 


4 


4.5 


52 


IGM Today: DBR 


151 




100 


760 


Solar system: 










during X-ray flare 


1.8 x 10 


-5 


3.5 x 10 4 


3.9 x 10 4 


from GDR 


1.0 x 10 


-5 


0.65 


12 


Quiet Sun 


1.7 x 10 


-5 


12 


18 


Recombo with CMB 


2.6 x 10 


-7 


8.0 


10 


Today from CMB 


92 




8.7 x 10" 3 


0.012 



IV, we give results for the time te 9 (M) (in years) for 
Zjv to reach equilibrium. Table IV goes on to give 
the binding energy of the least tightly bound electron, 
E B (M) for M > l(T 10 g and forM = lfT 21 g — for 
the seven settings. Recall from subsection II. B above 
that Eb = — Zjve 2 /6; 6 =max[rAr, ob/Zat]; 10 _10 g is the 
boundary, roughly, between the two regions. TE q (M) 
varies as M -1 / 3 for all settings. 

Other variations with ~ A/ +1 / 3 cancel out in 
the ratio Z^/Z: Fig. 4 does not have the kinks that 
are present in Fig. 3. The lack of any variation in 
E B (M) for M > 10" 10 should aid in devising SQN de- 
tection schemes, in the case of fragments from collisions 
or other production mechanisms giving a distribution of 
masses. The whole range would give the same signal. 
In those cases in which the charge focusing parameter 
f e of Eq. © satisfies f e >> 1, we can write a simple 
closed form for r Bq - Letting b =max[rjv, gb/Zn], writ- 
ing £ = 1(2) when b ^ cib/Zn {b — clb/Zn), and re- 
calling Eb = Z N ahc/((b), we have f e = AC,bE B / '(r N E e ). 
Consequently, it follows from Eq. ^ and the assumption 
that the integral in the capture rate, Eq. ([5]), is n e v e f e 
that we have 



Z N /Z N = nav e n e r N (hc/E e ) 



(9) 



Eq. ([9]) shows that the time to reach equilibrium depends 
only on the electron bath (or shower) and the SQN mass, 
but not the radiation; this is confirmed in the numerical 
results of Table IV and Fig. 4. 

We can go on to find a simple closed form equation for 
Eb, if b = r^v, in the approximations of either f e 3> 1 or 
f e < 1. Using Z + = Z- with Eqs. © and gives, for 
the first of these two cases 



F{E B ) = ir(v e n e r 2 N /(bE e ) 



(10) 



line width is from the sizes of f e which are, for cases 1, 3, 
and 5 respectively, 4.9 x 10 4 , 17, and 0.91. This gives at 
M = 1 g, r = 3.9 x 10~ 6 , 3.7 x 10~ 6 , and 2.2 x 10~ 6 years. 
Similarly, the two intergalactic settings (settings 2 and 
7), with DBR and CMB radiation, have f e = 9.9 x 10 3 
and 0.85 giving times of 33 and 20 years. Finally, COG 
today and universe just before recombination (settings 4 
and 5) have n e /v e = [40cm- 3 /60km/s; 287/300] giving 
f e = [44; 31] leading to r = [8.1 x 10~ 8 ;5.5 x 10" 8 ] This 
last case is "coincidental" as far as one can tell. 

In the calculations we first took the DBR and the 
galactic diffuse radiation (GDR) from Stecker and Sala- 
mon [23| combined with the graph of Henry [22j . How- 
ever, the result was binding energies only in the UV re- 
gions, so we switched to the reliable classic [28| (page 143) 
to generate our results. We said that the GDR around 
here (i.e. st the distance of the sun from the COG) could 
be extrapolated to the COG, which we took arbitrarily 
at 100 pc, by the r" 2 rule. 

Along with Fig. 4, in the second column of Table 



where F(Eb) is a function of a form that depends on the 
nature of the radiation distribution N y (here, thermal or 
power law). It is, for the case of f e ^> 1, the integral of 
the photon distribution function divided by Eb from Eb 
to infinity. 

Variation of Eb as a function of SQN mass M is small 
for small M, as well as absent for large, in the case of ther- 
mal radiation. This might be expected since the thermal 
spectral energy density cutoff is exponential (with Eb in 
the exponent). But it is significantly less rapidly varying 
for the diffuse radiation backgrounds where it is only a 
low power and hence Zn is more sensitive to the elec- 
tron bath or shower.. A similar separation of Z + = Z- 
into an equation with one side depending only on the 
radiation and the other depending only on the electrons 
can be made if f e <C 1. We do not provide estimates of 
Eb{M) for M in the transition region between small and 
large mass, since these should be most sensitive to more 
accurate wave functions of [|| . 

Some features of the results include the following: 
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In Fig. 3 we have used the parameterization of 
Eq. ||HJ| for the diffuse background radiation, de- 
rived from that plotted in [28( (intcrgalactic, in the 
area of the solar system same order of magnitude, 
and r~ 2 extrapolation to the COG). These are the 
only three non-thermal baths (or showers) consid- 
ered. 

The fact that Z N {M) behaves as M 1 / 3 , for the re- 
gion in which b = rjv, i.e. where M is large enough 
that the Bohr radius is inside the SQN, can be seen 
from Eq. 6. lib = rjy holds, all explicit dependence 
on M disappears from that equation and one just 
solves once for Eb = ahZ^/b, the binding energy 
of the most loosely bound electron. That the de- 
pendence of Zn on M is the same as that of b 
implies Z N ~ M 1 / 3 . 

Physically, one expects that the value of Zn should 
rise with M so as to make the binding energy of 
the most loosely bound electron independent of the 
radius of the SQN in the M region for which irb 2 = 
Ttr 2 N . The more interesting question is that of the 
negative slope of the little tails on the left in Fig. 3. 
(The answer is that, in that M region, Z_ grows 
as M 2 / 3 from the increase in r^r, while Z+, with 
b fixed at the appropriate Bohr radius, lacks that 
increase. This drives the solution to Eq. (J7J) to lower 
Zn to compensate.) 

The largest Zn values occur for solar system X- 
ray flares. Since these only last for minutes, one 
will need to consider carefully whether there is suf- 
ficient time to realize the large values and, even 
more importantly, sufficient time to exploit them 
for SQN detection. One sees, from Fig. 4, that the 
SQN mass must be over a ton before the time to 
reach equilibrium gets down to minutes. 

The results for SQNs at recombination are about 
the same as the results for the (quiet) sun because, 
while the temperature in the sun is higher than that 
in the CMB just after recombination, the electron 
density then was lower than that in today's solar 
wind near the Earth. 

Extragalactic SQNs have relatively high Zm values 
and ones near the solar system much lower ones 
because of the great difference in electron densities. 

Note that Zn would be limited by rapid vacuum 
breakdown to less than Eb(Zn) = —2m e . For 
Eb < -2m e , it is energetically favorable for the 
vacuum to create an electron with Eb = —2m e 
with one m e going to the mass of the electron 
and the other m e used to make a positron which 
goes off to infinity. However, Madsen [H| has re- 
cently studied, in some detail, pair creation impli- 
cations for SQNs. He points out that the rate of 
pair creation can be very slow, because the positron 



must tunnel through a Coulomb barrier. From his 
Eq. (10) one sees that the time for pair creation 
rises like M 1 / 3 and passes one second at SQN mass 
of about one gram. Thus the 2m e gammas will 
likely be at a higher energy for which the ioniza- 
tion rate equals the electron capture rate. For 
10 _20 g< M < 10 _18 g, pair production occurs at 
Z N > 137. For M < 10~ 20 g, there arc not enough 
electrons for it to occur at all. We discuss detec- 
tion of the 2m e line further in the last subsection 
of Section IV. Three of us, in a separate publica- 
tion, [39| have precisely computed the critical Z, as 
a function of SQN radius, by computing the Dirac 
Equation wave functions for E = — 2m e for two 
different charge distributions. 

• Our results do not take account of the energy cost 
of charge separation or of the stabilizing force of 
surface tension. The net effect of these two factors 
can be to fractionate large SQNs if surface tension 
is small or to disfavor small SQNs if it is large. 
The literature does not yet seem to contain reliable 
calculations of these effects, or the surface tension 
itself. In this paper we have considered the charge 
that would result for the full range of mass values, 
in spite of the fact that some values may turn out 
to decay into smaller SQNs rapidly or even to be 
precluded entirely. 

• The Thomas-Fcrmi approximation is valid where 
the electron wavelengths are smaller than the dis- 
tance, rjv, over which the quark electric field varies 
appreciably. Thus the transition with increasing 
r^v to the Thomas- Fermi regime might take place 
at an rjy such that the electron wave function be- 
comes appreciable in or near the quark nugget — 
around r^r ~ 10~ 10 cm. This is just one order of 
magnitude smaller than the Bohr radius of the least 
bound electron. References [H| and [39[ compute 
the transition to the Thomas-Fermi regime more 
precisely, the latter by explicitly solving the Dirac 
equation. 

IV. DISCUSSION 

There are many detection techniques to consider once 
the order of magnitude of the charges on SQNs, and 
the binding energies of their most loosely bound elec- 
trons, are known. One could look for emission lines 
from electrons being captured by SQNs. One could look 
for absorption lines and/or absorption edges in radia- 
tion coming toward us from behind an SQN population. 
One could try to detect the static charge along with the 
very small charge-to-mass ratio, of nearby SQNs passing 
through the solar system. Finally there is the possibility 
of indirect detection by means of some astrophysical ef- 
fect of a population of charged SQNs. Below we give a 
few examples of such lines of inquiry. The main problem 
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in looking for a signal will be the low SQN abundance 
given the number density limit pdm/M > tisqn(M). If 
we got lucky, there could be a large SQN shower from 
a relatively nearby and recent collision of the neutron 
stars in a binary system, but it would be difficult to ad- 
vise, with a straight face, anyone to spend real money on 
such a possibility. 

Particle Detection Techniques. One could consider space- 
based particle detection techniques, including ones that 
might be based on the Moon. The major capability 
needed is wide area coverage. Consider the event rate 

dN ev /dt = hsqnvsqnA (11) 

where n and v are the number density and speed of the 
SQNs and A is the effective detector area. Again, suppose 
all SQNs are of mass M and that usqn ~ Pdm/M- If 
we want to detect one SQN in the time t, we need (recall 
Pdm ~5x 10~ 25 g/cm 3 ) At/M > 10 17 assuming that 
vsqn ~ 250km/s, the galactic virial velocity. If a square 
kilometer could be instrumented, nanogram SQNs might 
be detected at rates up to 100/s, and one-gram SQNs 
once a year. Something to consider is whether SQNs 
entering the atmosphere would produce a characteristic 
scintillation signal detectable by ground or space-based 
telescopes [40| . 

Absorption and Emission Lines and Edges. A second 
approach is to search settings in which there is identifi- 
able absorption or emission. One example would be to 
look for settings in which there is sufficient high energy 
radiation to bring SQNs to a high enough degree of ion- 
ization that pair creation ensues if another electron is 
ejected. At zero temperature, we expect that to happen 
when Eb, the binding energy of the least bound electron, 
is 2m e c 2 . This follows simply from energy conservation 
(E B + m(e+) + m(e") = 0). See Madsen [H for dis- 
cussion of other properties of low temperature SQN pair 
production 

At such a point, there should be an emission line at 
E = 2m e c 2 from electron capture into those states with 
energies near Eb = —2m e . This would result from a 
strong enough photon distribution to preclude lower Eb 
while, at the same time, vacuum pair creation prevented 
a larger one. The upshot would be that any electron 
capture would produce a 2m e photon, a very distinctive 
signal. However, as noted above, Madsen shows that 
the time required to produce a minimum energy pair 
rises exponentially with M 1 / 3 for SQNs, passing 1 sec- 
ond around a gram and shifting the 2m e signal toward 
higher energy. We emphasize that careful calculations 
are needed to determine whether that signal would be 
distinguishable against the strong backgrounds in an ex- 
plosive event. Once Eb is at a high enough energy to 
permit rapid pair production, any further electron ejec- 
tion results in an e + e~ pair so that there would also 
be gammas with E^ = m e c 2 from e + — e~ annihilation, 
where the positron annihilation could be inside the SQN 



or outside the SQN with an ambient electron. This signal 
should occur as the SQN, in the radiation and electron 
baths, oscillates about the equilibrium value of Zn- It 
requires an SQN mass above about 10 -20 g so that the 
total number of electrons is sufficient to reach that point. 

In addition to gammas from the positrons annihilat- 
ing with ambient electrons and the 2mc 2 emission there 
might be an accompanying absorption edge at 2mc 2 since 
photons of energy over that value could liberate electrons 
from the SQN, but less energetic ones could not. Our 
preliminary calculations summarized above indicate the 
diffuse radiation background is not able to ionize SQNs to 
that degree. Other places to look would be at explosive 
events, including supernovae (of various types), gamma 
ray bursts, and neutron star bursts and superbursts. 

It is possible that GLAST data, when available, could 
contain some evidence of this triad of signals: an emission 
line at E = m e c 2 , an emission line at E = 2m e c 2 , and an 
absorption edge starting from E = 2m e c 2 . The emission 
line at a single photon energy of 2m e c 2 (resulting from 
7's keeping Eb > 2m e c 2 and vacuum pair creation keep- 
ing Eb < 2m e c 2 ) would be a dramatic SQN signal. It 
may well be that the background would swamp any sig- 
nal, but this three-pronged test could be such a definitive 
indicator that it is important to assess its feasibility. We 
know of no other source of 2m e c 2 emission lines. A place 
to look for an analogous signal, without the emission line, 
would be in the solar system during a solar X-ray flare. 
Then we would expect an X-ray signal of about 50 KeV 
and some absorption of X-rays over that (see Table IV 
and Fig. 4). Additionally, we note that a strong X-ray 
signal from electron-positron annihilation strongly con- 
centrated at the center of our galaxy has been observed 
most recently with the INTEGRAL satellite. See, for ex- 
ample, Yuksel [4l[ for a brief review. We have, however, 
not seen any reports of 1.02 MeV photons, so it would be 
premature for SQNs to clamor for entry into Yuksel's list 
of about two dozen "exotic" models that might account 
for the positrons. 

A second place pair production is thought to have oc- 
curred was in pair instability supernovae [42|. At least 
one recent supernova, 2006GY, discovered by Quimby 
et al. [431, and studied by Ofek et al. and Smith et 
al. [!U showed strong positron production. It was the 
most luminous supernova recorded and had a light curve 
that lasted significantly longer than the typical super- 
nova. Leahy and Ouyed [45| have explained its features 
on the basis of a two step process, the first an ordinary 
type II supernova implosion followed by a transition from 
a neutron star to an SQS. A second explanation might 
conceivably be possible using results along the lines of 
Fig. 4. A cloud of gram range SQNs from an event with 
persistent production of a lot of MeV range 7s would gen- 
erate positrons. It could also generate the 2mc 2 j signal 
discussed above and not expected from any other obvi- 
ous, at least to us, processes. This could be the proverbial 
"smoking gun" for the presence of SQNs. Pair produc- 
tion as a cooling mechanism for hot SQSs was pointed 
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out by Usov [3l| and is an active area of study [32| . 

Other Settings. In addition to the 7 settings treated here, 
and the explosive case cited, there are many other places 
where electrons can be stripped from SQNs. One exam- 
ple is SQNs captured into the Earths geomagnetic field. 
Paulucci, Horvath and Medina- Tanco [46| have consid- 
ered this setting with trapped electrons both ionizing and 
de-ionizing the SQN. Another interesting setting would 
be an SQN passing through the sun. There, collisions 
with photons and protons would ionize: protons are much 
more effective than electrons because of the larger mass: 
collisions with an SQN with the galactic virial velocity 
gives a fraction of an eV for electrons compared to a KeV 
for protons. The balance between ionization by photons 
and protons and electron capture by collisions with elec- 
trons should give a relatively high ionization, especially 
if the still greater energies possible from collisions with 
helium are taken into account. The question would be 
one of the distance from which such a charge could be 
detected. Passage of an SQN through a giant molecular 
cloud might be of interest. 

Another important setting would be passage of an SQN 
through metal, rock or other common materials. Esti- 
mating the charge produced by passage through rock, 
for example, would help determine whether detection in 
a deep underground laboratory would be feasible as well 
as the sort of signal that other detectors might give if a 
highly charged, massive body transited them. 

In summary, we have computed, for seven different set- 
tings, the charge that would accumulate on an SQN pass- 
ing through a photon bath (or shower) with ambient elec- 
trons. In Fig. 3 we give the resulting equilibrium charges 
as a function of SQN mass. The largest charge among 
the settings considered for an SQN passing through the 
solar system during a solar X-ray flare. In Fig. 4 we 
give an estimate of the times required to reach the equi- 
librium values. In Table IV we give the corresponding 
equilibrium binding energies of the least bound electron. 
We give these results for SQNs of masses ranging from 
10~ 20 to 10 +30 g. 

V. ARE THEY THERE? 

It is now three dozen years since Bodmer's paper Q 
and two dozen since Witten's [l|, but there has been no 
experimental or observational sign of SQM - in spite of 
a reasonably vigorous program of searching. One has to 
ask the questions both of whether the game is worth the 
candle and of what is the score anyway? 

First the score. There are no astronomical sightings 
of strange quark stars (SQS 's). As noted, two reports of 
strange quark pulsars have gone away [l?]] . Addition- 
ally, there are two sets of observations that appear to pre- 
clude the neutron stars observed being SQSs. First, there 
is the glitch problem. Pulsars, from time to time, sud- 
denly change rotation frequencies. This change is most 
easily interpreted as the result of starquakes in the less 



dense, more familiar matter on top of the neutron core. It 
is difficult to explain these glitches with SQSs since there 
is a limit to the mass of a familiar matter crust that can 
rest on a SQM core [47| (although one SQM model ap- 
pears to do so (HI). Then there are the explosions from 
time to time on soft gamma ray repeaters (SGRs). These 
are believed to be neutron stars: some bursts are quite 
large and are, of course, called superbursts. The X-rays 
from superbursts show a structure of quasi-periodic os- 
cillations (QPOs) during the course of the burst. QPO 
frequencies identified during a couple of relatively recent 
superbursts correspond to those expected from a thicker 
crust [49l | than neutron drip permits an SQM core to 
maintain [2^ . 

Moving on to terrestrial (and lunar) searches, observa- 
tions and experiments were noted in Section I and sum- 
marized in Table I. A more thorough review is given in 
[ll[, references therein and other papers from that con- 
ference. 

Last but not least is theory. Recognizing that Cooper 
pairing should add tens of MeV per baryon to SQN bind- 
ing energies 0] greatly enhanced the argument for SQM 
existence. On the other hand, while observations have 
been casting doubt on the existence of SQSs (and frag- 
ments of SQS - SQS collisions), theory has been casting 
doubt on primordial SQN production. The current situ- 
ation is summarized by Boyanovsky et al. fic| along the 
following lines: 

There is an estimate that a primordial SQN must have 
mass M greater than 10 20 g in order for neutrino cooling, 
which scales as M, to be greater than evaporation cooling 
which scales like M 2 / 3 [? ]. Furthermore, if the (first 
order) phase change to SQM occurs, as expected, when 
the strong interactions become strong at about 200 MeV, 
one would expect a mass of this same order - give or take 
a factor of 10 3 or so. This is about the mass of a large 
comet. Assuming the solar system formed from a light 
year of gas, at a density of about 5 x 10 _25 g/cm 3 , the 
galactic dark matter density in our neighborhood, there 
would, on the average, be one SQN every 10 15 cm (60 
AU). There could be on the order of a billion SQNs in the 
solar system, some in the Oort cloud, some in the Kuiper 
belt, and some in the interiors of the sun and the planets, 
but not many in nearby orbit about the sun since orbit 
perturbations by Jupiter and other planets tend to drive 
comet-sized objects out of the inner solar system. There 
should also be SQNs passing through the solar system if 
they exist and if they constitute the galactic dark matter. 

These would seem clear indications for where to con- 
centrate future searches for primordial SQNs, but there 
appears to be a fly in the ointment: again the neutri- 
nos. By diffusion, they wipe out inhomogencitics over 
distances larger than the neutrino mean free path in the 
time it takes for hadron bubble nucleation and quark- 
gluon plasma reheating. Returning neutrinos cancel neu- 
trino emission. 

Are they there? Is the game worth the candle? Finch 
says in [1J| "If it turns out that SQM is stable, the impli- 
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cations would be potentially tremendous not only for the 
resultant direct and indirect understanding of the strong 
interactions but also for the practical applications rang- 
ing from new materials..." In brief, the answer to "are 
they there?" seems to be that it's a long shot, but to 
some worth trying. 

Our hope is that finding the charge on SQNs in space 
will help devise space tests that either detect SQNs or 
narrow the mass windows permitted. 

Acknowledgments 

VLT very much appreciates a number of very help- 
ful conversations with Demos Kazanas, one with Floyd 



Stecker on reference [23| and one with M. Alford effects of 
Cooper pairing. Dr. Alford kindly provided a number of 
very helpful comments on the first version of this paper. 
DCR and VLT are also grateful to Jonathan Gardner 
for calling SN-2006gy to our attention, as well as C. Kil- 
bourne for talking WHIM with us. DAD was supported 
in part by the U.S. Department of Energy under Grant 
No. DE-FG03-93ER40757. WWR was supported in part 
by the National Science Foundation under Grant PHY- 
0555544. 



[1] E. Witten, Phys. Rev. D 30, 279 (1984). 
[2] A. de Rujula, and S. Glashow, Nature (London) 312, 734 [25 
(1984). [26 
[3] In preparation. 
[4] A. R. Bodmer, Phys. Rev. D 4, 1001 (1971). [27 
[5] E. Farhi and R. Jaffe, Phys. Rev. D 30, 2379 (1984). [28 
[6] A. Chodos el al, Phys. Rev. D 9, 3471 (1974). 
[7] M. Alford, K. Rajagopal, and F. Wilczek, Phys. Lett. [29 
422B, 247 (1998). [30 
[8] M. G. Alford, A. Schmitt, K. Rajagopal, and T. 

Schafer, Color superconductivity in dense quark matter [31 
larXiv:0709.4635l [hep-ph] (2007). 
[9] J. Madsen, Phys. Rev. Lett. 87, 172003 (2001). [32 
[10] D. Boyanovsky, H. J. de Vega, Ann. Rev. Nucl. Part. Sci. 

56, 441 (2006) ; [a^v:hep-ph/060"2002] v2 27 Jul. 2006. 
[11] E. Finch, "Strangelets: Who is looking, and how?" at [33 
SQM2006: Conference on Strangeness in Quark Mat- [34 
ter, Los Angeles J.Phys. G32:S251-S258,2006. e-Print: [35 
|nucl-ex/0605010| 
[12] J. Sandweiss, J. Phys. G 30, 51 (S) (2004). [36 
[13] P. Guinnessy, Physics Today 60 no. 5, 30 (2007). 
[14] C. Spiering, Proceedings of ICRC2001 1242, (2001). [37 
[15] P. B. Price, Phys. Rev. D38, 3813 (1988). 
[16] J. J. Drake et al, ApJ. 572, 1996 (2002); P. G. Slane et [38 
al, ApJ. L45, 571 (2002). [39 
[17] F. M. Walter and J. Lattimer, ApJ. 576, L145 (2002). 
[18] W. B. Banerdt, T. Chui, E. T. Herrin, D. C. Rosenbaum, [40 
and V. L. Teplitz, in IDM 2004: 5th International Work- [41 
shop, on the Identification of Dark Matter, Edinburgh, [42 
edited by V. Kudryavtsev (World Scientific, New Jersey, [43 
2005) p. 581. [44 
[19] E. T. Herrin, D. C. Rosenbaum, and V. L. Teplitz, Phys. 

Rev. D 73, 043511 (2006). [45 
[20] J. Horvath, Astrophysics and Space Science DO 10 [46 
(1997). 

[21] A. Zhitnitsky, Phys. Rev. D74:043515 (2006). [47 
[22] R. C. Henry, ApJ. 516, L49 (1999). [ 
[23] F. W. Stecker and M. H. Salamon in Gamma 2001, edited 

by S. Ritz, N. Gehrels and C. Schrader (AIP, New York [49] 

2001) p.432. 

[24] M. G. Alford et al. , Phys. Rev. D73, 114016, 2006 (arXiv 



hep-ph/06 04134). 

M Oertel and M. Urban. larXiv:0801.2313"V 2 [nucl-th]. 
J. Binney and S. Tremaine, Galactic Dynamics (Prince- 
ton University Press, Princeton, 1988). 
V. Safronov, Ann. d'Astrophysique 23, 979 (1960). 

E. W. Kolb and M.S. Turner, The Early Universe 
(Addison- Wesley, Massachusetts, 1990). 

C. Alcock, E. Fahri, and A. Olinto, ApJ. 310, 261 (1986). 
L. I. Schiff, Quantum Mechanics 3rd edition ( McGraw- 
Hill, New York, 1968) p. 427. 

V. V. Usov, Phys. Rev. Lett. 80, 230 (1998) and 81, 4775 
(1998). 

See, e.g., T, Harko and K.S. Cheng, ApJ. 643, 318 
(2006); S. P, Kim, hep-ph 0801.3300 and references 
therein. 

J. Schwinger, Ph ys. Rev. 82, 664 (1951). 

F. Paerels et al. , larXiv:0801 . 10641 [astro-ph] 7 Jan 2008. 

D. Scott et al, in Allen's Astrophysical Quantities, edited 
by A. N. Cox (Springer, New York, 2000). 

B. W. Carroll D. A. Ostlie, Introduction to Modern As- 
trophysics, (Addison- Wesley, Massachusetts, 1996). 

E. Komatso et al. arXiv:astro-ph/0803.0547 submitted to 
ApJ. 

J. Madsen, Phys. Rev. Lett. 100, 151102 (2008). 

D. Dicus, W. Repko, and V. Teplitz arXiv 0808.1739 
(2008) (hep-ph), Phys. Rev. (to be published). 

R. Schwitters, private communication. 

H. Yuksel, [asFro-ph/0609139| 

W. Fowler and F. Hoyle ApJ. 9, 201 (1964). 

R. M. Quimby et al, ApJ. 668, L99-L102 (2007). 

E. Ofek et al, ApJL, 659, 13 (2007) and N. Smith et al, 
666, 111 (2007). 

D. Leahy and R Ouyed, astro-ph 0708.1787 (2008). 

L. Paulucci, J. Horvath and G. Mdina-Ranco, Phys. Rev. 

D77: 043003 (28). 

F. Weber et al, arXiv:astro-ph/0705.2708. 

P. Jaikumar, D. Reddy ,and A. W. Steiner, Phys. Rev. 
Lett 96, 041101 (2006). 

A. L Watts and S. Reddy, MNRAS L, 63, 379 (2007). 



